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Abstract. The microcanonical dynamics of an ensemble of random magnetic dipoles 
in a needle has been investigated. Due to the presence of a constant of motion in the 
1-D "dimensional" phase transition in the quasi one dimensional case has been 

found separating a paramagnetic chaotic phase from a ferromagnetic regular one. In 
particular, a simple criterium for the transition has been formulated and an intensive 
critical parameter found. Numerical simulations support our understanding of this 
complex phenomenon. 
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1. Introduction 

A truly comprehensive understanding of magnetism at the nanoscale is still lacking and 
has important consequences in the technology of memory and information processing 
devices. 

Many unsolved problems about magnetic properties of diluted spin systems 
attracted recently great attention. Among the open problems there is the emergence of 
ferromagnetism in doped diluted systems[I], where the Curie temperatures can be as 
high as 300 K, and a deep theoretical understanding of the magnetic properties of dilute 
dipole systems (spin glass transition, ferromagnetic and anti- ferromagnetic transitions). 

Here we will concentrate on randomly arranged dilute classical dipoles, which are 
called dipole glasses. Many results in literature, sometimes controversial, exist on such 
kind of systems. Magnetic properties of dipole-dipole interacting spins are particularly 
difficult to study due to many factors: long range nature of the interaction, anisotropy 
and frustration. Long range and anisotropy can induce ergodicity breaking [2] in a 
system. Breaking of ergodicity, a concept introduced by Palmer [2], and recently found 
explicitly [H H] in a class of long-ranged anisotropic spin systems, is a key word to 
understand phase transitions too, even if it should not be confused with breaking of 
symmetry [5]. Speaking loosely, few constants of motion, such as the energy, or the 
angular momentum, in a particular geometry, produce a separation of the allowable 
phase space in two or more subspace over which the motion is constrained. In Ref . [3] the 
energy at which the separation occurs has been calculated explicitly for an anisotropic 
1-D classical Heisenberg systems. In that case both the anisotropy and the long ranged 
nature of the inter-spin interaction, are essential ingredients in order to have breaking 
of ergodicity [7]. On the other hand, frustration, that is the impossibility to attain a 
global minimal energy minimizing locally the interactions, induces a dependence of the 
ferromagnetic and anti- ferromagnetic properties on the lattice geometry [8]. 

Other results concerning the so-called Ising dipole glass can also be found in 
literature, where Ising simply means uni-axial. To quote but a few: spin glass transition 
for high concentration, using Monte Carlo simulation[9l [10], mean field spin glass 
transition at low concentration depending on the lattice geometry [11]. absence of spin 
glass transition for low concentration using Wang-Landau Monte Carlo simulations [12] 
or the recent spin glass transition at non zero temperature from extensive numerical 
simulation [T3]. 

In this paper we will focus our analysis on a dipole glass of freely rotating classical 
dipoles. First of all the dipole glass is a typical example of very frustrated system 
[TH [T5l IT6] . so that different ground state configurations can exist depending on the 
geometry and the spin concentration. Results in the canonical ensemble typically 
consider a mean field approach, and it is common lore that the random positions of 
the dipoles induce magnetic field fiuctuations. These fiuctuations do not vanish at 
T — 7- 0, unlike thermal fiuctuations, and tend to suppress magnetic order even at 
T = 0[T5l [16]. So, magnetic order, is expected to happen only for high impurity 
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concentration (and small temperature) [TBI [13 EH]- Mean field theories consider only 
the equilibrium properties and do not take into account the time needed to reach the 
equilibrium situation and finite size effects. On the other hand the question of how long 
a metastable state can last is a major issue in determining the magnetic properties of a 
system. 

In this paper we study the microcanonical dynamics, reserving the study of the 
infiuence of a thermal bath for further investigations. We analyze the microcanonical 
dynamics of dipoles put at the vertexes of a cubic lattice (so that their relative distance 
cannot be smaller than the lattice size), only on the basis of the Landau-Lifshitz-Gilbert 
equations of motion. 3-D dipole-dipole interacting systems can be realized quite easily 
in laboratory, for instance doping a non magnetic media with paramagnetic ions, weakly 
interacting with the lattice and with an inter-dipole distance sufficiently large in order to 
neglect Heisenberg interaction ( low concentration). The choice of studying a random 
glass instead of a system composed of dipoles regularly arranged in some lattice is 
twofold: from one hand is it relatively easy to dope a system putting some paramagnetic 
doping ions in a random way inside any non-magnetic media. On the other hand, a 3-D 
cubic lattice with a full concentration of dopant ions 6 = 1, even if thin, does not have a 
ferromagnetic ground state [19J, so that another type of transition should be considered 
(paramagnetic / antiferromagnetic) . 

Anticipating some of the results, we have found that taking into account a typical 
experimental situation with needle-shaped sample, at low dipole concentration a further 
constant of motion appears that induce a kind of "phase" transition related to invariant 
tori, which separate the allowable phase space in many disconnected regions. This 
result seems to indicate that at very low concentration a system of random dipoles 
in a needle resemble a one dimensional arrangement of dipoles, and thus can have 
ferromagnetic behavior. In this particular case, the ergodicity breaking is not due to an 
increase of energy, but to an increase of perturbation, which means the tendency to a 
transformation from a needle shape (quasi 1-D system) to a cubic shape (3-D shape). 
In a sense, these results are more akin to the standard perturbation theory in classical 
dynamical systems |20] |2T] . re- interpreted in the light of phase transitions induced by 
demagnetization times |22] . 

In the future we are going to study the same system in contact with a thermal 
bath. In this case the presence of the ergodicity breaking found in [3l H] should infiuence 
the demagnetization times. In the microcanonical case, the presence of this ergodicity 
breaking is hidden by the quasi-integrability of motion. 

2. The Classical Model and the Perturbative Approach 

Let us consider a system of classical dipoles /Ij randomly put at the nodes of a 3-D 
gridded box R x R x L, with L ^ R, and low concentration 5 -C 1, as indicated in 
Fig.m 

From the physical point of view it represents a dilute system of paramagnetic ions in a 
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Figure 1. Needle geometry. The classical dipoles are put in a random way on the 
vertexes of a cubic lattice of size a. R and L are given in units of the lattice size a. 



non magnetic bulk, with a concentration 6 = N/Ng where Ng = R^L is the number of 
allowable sites in the 3-D lattice. As explained above, such a system can be realized 
in laboratory, doping a non-magnetic system having a cubic lattice with paramagnetic 
impurities (e.g. doped Ti02 and other [1]). If the dipoles weakly interact with the lattice 
and if their average distance is much greater than the Bohr radius, we can simply 
neglect the Heisenberg (exchange) interaction and represent their mutual interaction 
and dynamics with a pure dipole-dipole interaction energy: 

,, ,,2 1 

i=l j>i \' %j\ 

Here Si is the i-th dimensionless spin vector 

S^■S,= 1, (2) 

/i is the magnetic moment of the paramagnetic doping ions and r^j is the distance 
between the i-th and the j-th spin in units of the lattice spacing a. 

The dynamics is described by the Landau-Lifshitz-Gilbert equations of motion: 
d ^ ^ 6E 

^/.. = 7/^.x— , (3) 

where flk = fiSk and 7 is the gyromagnetic ratio. 
They can be rewritten in the dimensionless form. 
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where the following dimensionless quantities have been introduced: 
Eq = E 

T = ut, With 00 = - — -. 

The system of equations (j4]) conserves the energy ([1]) and the squared modulii of 
the spins ([2]). 

The diluted quasi 1-D system can be magnetized with a strong magnetic field 
directed along L, the longest axis (2;-axis). The questions we would hke to answer is 
the following: What is the dependence of the average demagnetization time and its 
fluctuations on the system parameters? 

The relevant parameters to take into account are the concentration S of 
paramagnetic ions and the aspect ratio e = R/L. In principle, due to the long-ranged 
nature of the dipole interaction, one could ask whether there are effects dependent on 
both the system size and the number of doping spins A^, even if in quasi 1-D systems, 
the dipole interaction can be treated as a short range interaction. 

From the point of view of the equations of motion (jlj), if the dipoles are lying 
along a straight line (i? = ^ e = 0), there is a further constant of motion, i.e. 
Mz = J2k ^k- Therefore, for a 1-D system, the answer to the first question above 

is very simple : a state with any initial magnetization Mz{0) 7^ will keep the initial 
magnetization forever. The natural question thus becomes: what happens for e 7^ 0? 
Will a magnetized state demagnetize and how much time it takes to do that? 

The classical dynamical picture can be simplified adopting a perturbative approach, 
namely approximating the unit versor between two spins as follows: 

cos 9ijZ + sin 9ij (cos (pijx + sin (pijy) ~ 5 + (eA^) (cos (pijX + sin 4>ijy) (6) 



where x, y, z are the unit versors, (pij are the azimuthal angles with respect the pr- 
axis and 9ij are the polar angles. In the last equation we approximate cos 9ij ~ 1 and 
sin 9ij ~ R/ (d), where, for dilute dipoles in a needle goemetry, (d) ~ L/N is the average 
distance among spins. The energy ([1]), to first order in eN, becomes: Eq = Hq + eNV, 
where Ho is the energy part that conserve M^, and V is the perturbation, 

N 1 



TT ^ \ ^ \ ^ QX QX : Oy Oy O QZ QZ 

-"0 — 7:2^2^] 13 iJi >Jj -r 'Jj '^iJi'Jj: 



2 i=l j^i Vij\ 
N 



(7) 



y = -3EE 



1 j^i 



|3 



COS (pijS-Sj + sin (pijS-Sj 



The equations of motion for the macroscopic variables, M^^y^z can be written as, 
^7-^ = 3e E r— nr'S'f ('5'^ cos (pik - SI sin 0^^.) 

^ = ^ E e:^ {StSl + eN [SlS^smct^^u + {SIS^ - SlSt) cos0,,]} 
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Figure 2. Three different trajectories of the total magnetization at the plane Mz = 0, 
for i? = 4, L = 4000, 6 — 10"'^, N — 64, in the integrable case. Initially spins are 
chosen with random components on the unit sphere. 



dr 



- If E T^iSlSt + eN [S^Sf cos 0.. + {StSf - S^S^ 



sm 



and, in particular, for e 

dM, 



we have : 



dr 



(9) 



dM^ 
~d^ 



k 



having defined, the average "local" frequencies: 

1 



3E 



rik\ 



(10) 



These equations describe a kind of rotation in the plane perpendicular to the z- 
magnetization (which is a constant of motion). Therefore one could expect that for 
eN <^ 1 a rotational-like motion about the 2;-axis persists, while remains a quasi 
constant of motion. This is what can be observed for instance by a direct inspection 
of the trajectories of the macroscopic vector M, in the plane x, y, see Fig. O where 
few selected trajectories has been iterated in time, for e = 10~^ and N = 64. Quite 



Dynamics of random dipoles : chaos vs ferromagnetism 



7 



-0.4 



Paramagnetic 



Ferromagnetic 




-0.2 





AM 



0.2 



0.4 




Figure 3. Data in this Figure refer to systems with iV = 64 spins and a concentration 
S = 10^'^. The time behavior of the magnetization is shown, for different initial 
conditions, in the overcritical case e = 0.125 (a) (L = 160, i? = 20) and in the 
undercritical case e = lO^'^ (b) {L — 4000, i? = 4). In c) and d) the probabihty 
distribution functions for the fluctuations AM^ = i{M^) — (M^)^)-^/^ around the 
equilibrium value is shown for the data given respectively in a) and b) . 



naturally, on increasing the perturbation strength e, one could expect that the invariant 
tori Mz = const will be broken, and, eventually, a stochastic motion of the macroscopic 
variable will appear. In the next Section we will study the survival of invariant tori, 
under the dimensional perturbation eN > 0. 

3. The chaotic— paramagnetic and the integrable— ferromagnetic phases 

The dynamical behavior of the system can be characterized by a "regular region" eA^ < 1 
in which the magnetization Mz{t) is bounded in a small interval SM^, while, for eA^ > 1, 
Mz{t) quickly decays to zero. To be more precise, the transition across {eN)cr — 1 is 
smooth, namely there is a region of eN values in which the initial magnetization decay 
to some non zero constant when the time r becomes large. 

The critical value of the perturbation strength (eiV)cr — 1 can be obtained with 
the following hand- waving argument. Let us divide the 3-D box mn = L/R = l/e 
small cubic boxes of side R. If the impurities concentration 6 is sufficiently small in 
order to have only one spin inside each i?-side box then the system is approximately one 
dimensional and Mz can be considered an approximate constant of motion. Otherwise, 
for large 6, the system behaves like a 3-D system and can spread everywhere. In 
other words, in order to have less than one spin in each block one should have 
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N/n < 1, or eA^ < {eN)^r = 1- 

Moreover the study of the dynamics done in the previous Section suggests to take 
as a small parameter eN and to look for ferromagnetism when eN < 1. Note that this 
choice is also appropriate from the thermodynamic point of view since eA^ = RN/L is 
an intensive parameter in the large limit L — oo, N — )■ oo, N/L = const, with 
R fixed. 

An example is shown in Fig. |3l where the dynamics of magnetization has been 
plot in the overcritical case (eA^ > 1 Fig. [3|l) and in the undercritical one (eA^ < 1 
Fig. ISh"). Different trajectories, corresponding to different initial conditions Mz{0) have 
been shown in different colors. As one can see, in the "paramagnetic" phase (eA^ > 1) 
the magnetization first decays to zero and then it fluctuates randomly around zero. 
On the contrary, in the "ferromagnetic" phase (eA^ < 1), it shows a periodic behavior 
around the initial conditions. 

This behavior is quite typical in the study of dynamical systems, where the increase 
of a suitable perturbative parameter is related to the breaking of invariant tori and to 
the emergence of chaotic motion [201 121] . 

It is also remarkable to study the fluctuations around the asymptotic behavior: in 
the undercritical case (Fig. |3]i) fluctuations are much smaller than in the overcritical 
case (Fig. roughly 10 times for this case, as can be seen comparing the width of the 
probability distribution functions in Fig. |3]c) and d). 

The large fluctuations around the average values and in order to fit a possible 
experimental situation, suggest averaging over disorder, namely an ensemble of samples 
with different random configurations, initially magnetized along the 2;-axis. 

The results for the ensemble average (M^) are shown in FigJUi, where the different 
behavior in the two "phases" eA^ < 1 and eA^ > 1 are reflected in an average 
magnetization not decreasing or decreasing to zero. Indeed, the average magnetization 
in the undercritical regime reaches some equilibrium value different from zero after some 
initial decay, while in the overcritical regime it goes to zero in an algebraic way. 

Ensemble fluctuations at the equilibrium are independent of eA^ in the paramagnetic 
phase while in the ferromagnetic one are typically smaller and increasing as \feN . They 
are presented in Fig. HJd), where 

AM,^^ = lim AMJr) 

has been shown as a function of eA^. On the vertical axis we renormalize the asymptotic 
values by \J1^N), to take into account fluctuations due to variation of the number of spins 
A^. Each set of points on the plot corresponds to an ensemble of magnetized needles, 
with the same concentration 6 and number of spins A^ (paramagnetic ions) and different 
aspect ratio e, or same concentration and aspect ratio and different number of spins. It 
is quite remarkable that the critical value eA^ ~ 1 , is well fitted by all different series, 
suggesting eA^ as a good scaling parameter for the macroscopic behavior. 

Both the independence of the perturbation strength in the paramagnetic phase and 
the square root dependence on eA^ in the ferromagnetic phase can be understood on 
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Figure 4. a) Time beiiavior of tlie average magnetization for different values of the 
aspect ratio e = 1.25 x 10"*, 2.9 x 10"'', 9.8 x lO""*, 4.5 x lO'^, 2.5 x IQ-^, IQ-^ (from 
the upper to the lower), fixed number of spins N — 220 and fixed concentration 
6 = 10^'^. The average is taken over 100 different random configurations. Initially 
we choose S'f(O) = 1, i = 1, . . . , A^. b) Dependence of the equilibrium value of 
fluctuations as a function of eN. Dashed vertical line indicates the critical value 
eN = 1. Red line indicates the dependence V Ne. Different symbols stand for : 
TV = 40, 6 = 0.5 10-3 (circles), = 80, S ^ 10^^ (squares), N = 220, 6 = 10"^ 
(triangles down), e — 0.1, S — 10^^ (triangles up). 





the basis of classical dynamical theory. Breaking invariant tori with a perturbation 
strength k corresponds to create stochastic layers between invariant tori whose size is 
proportional to -\/A;[201 |2T]. On the other side when the system is completely chaotic, 
since the variable is bounded, it can only occupy all the allowable stochastic region, 
and a further increasing of perturbation strength can not modify this size. 

Finally we point out that around eA^ = 1 we can expect a transition from a 
ferromagnetic ground state to an antiferromagnetic ground state. Indeed for eA^ ^ 1 
the system is close to a ID arrangement of dipoles, so that the ground state will be 
ferromagnetic as pointed out in the previous Section. On the other side for eA^ ^ 1 
the system is close to a 3D arrangement of dipoles. In this case, for a simple cubic 
lattice, it can be shown [19j that the ground state is antiferromagnetic. Some numerical 
simulations we did confirm this conjecture, but work is still in progress and will be 
presented in a future publication. 

4. Demagnetization time 

As we have seen in the previous Section, the system dynamics can be described by the 
parameter, eA^, characterizing two different dynamical phases, and describing how much 
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Figure 5. Average Magnetization at time T as a function of cN. For the same eN we 
show the average magnetization after different simulation times, T. Data refer to the 
case 6 — 10^^. Different number of dipoles are shown: N = 220 (circles for T = 10'^ 
and lozanges for T = 10'^) and iV = 80 (squares for T = 10^, left trangles for T ^ 10^ 
and crosses for T = 10^). 



one-dimensional a system is. In this Section we will show that eA^ is also a good (and 
intensive) scaling parameter for the macroscopic properties of the quasi 1-D system. 

In order to prove numerically such argument we need to find two physical 
observables that can describe the paramagnetic and the ferromagnetic phase and to 
study their dependence on the parameter eN. To this end, let us introduce, on the 
paramagnetic side, the demagnetization time ri/2, defined as the time at which the 
average magnetization decay to one half of its initial value: 

(M,(ri/2)) = (1/2)(M,(0)). 

In the same way, on the ferromagnetic side, we introduce the "remnant magnetization" 
Mr as the magnetization left when r — ?■ cxd, 

Mr = Jim (M^(r)). 

Let us stress that both quantities are physically sound, in the sense that they are directly 
and easily measurable. 

While it is clear that, even if both quantities can be defined only in their respective 
phases, they can give useful information when extended to the other phases. For 
instance, in the paramagnetic phase the average magnetization will depend strongly 
on the simulation time (T), while in the ferromagnetic phase the demagnetization time 
is typically infinite|j| The dependence of the average magnetization at fixed time T can 

I In this case we underestimate its value putting the maximum dimensionless simulation time, which 
is 104. 
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Figure 6. Left axis : remnant magnetization Mr vs the parameter eN indicated 
witli open symbols. Rigth axis : inverse demagnetization time Ti/2 rescaled by the 
concentration 5, vs the parameter eN indicated with full symbols. List of symbols: 
circles (e = 0.1,5 = IQ-^), lozanges {N = 80,5 = 0.01), crosses {N = 220,6 = 10'^), 
squares {N = 100, i? = 4), asterisks (N = 500,7? = 20), left triangles (e = 0.01,5 = 
0.01). Initially we choose (0) = 1, i = 1, . . . , iV. In this Figure an ensemble of 100 
different configurations has been considered. Each member of the ensemble has been 
integrated for 10"^ dimensionless time units. 



give important information on the ferromagnetic-paramagnetic transition. Indeed we 
can expect a weak dependence of the average magnetization on the simulation time T, 
in the ferromagnetic region eA^ < 1, since the presence of quasi-constant of motions 
freezes the magnetization, while, in the paramagnetic phase, the average magnetization 
at time T, < Mz{T) >, goes to zero as the time grows. This fact is clearly shown in 
Fig. ([5]) , where the average magnetization < (T) > is plotted vs eN for different times 
T. As we can see from Fig. ([5]), as we increase the time T the average magnetization 
remains almost constant in the ferromagnetic side (eA^ < 1), while it goes to zero in the 
paramagnetic side (eA^ > 1), thus demonstrating a clear signature of the dimensional 
transition discussed above. 

Since eA^ = 6R^, one can essentially consider different ways to approach the critical 
point eA^ = 6R^ ~ 1, keeping fixed one of the four quantities e, A^, 6, R and varying 
correspondingly the others. This is exactly what we did in Fig. [6l where we show, 
on the same plot, the remnant magnetization (open symbols and vertical axis to 
the left), and the inverse demagnetization time, Ti/26, rescaled by the concentration 
S ( vertical axis to the right, full symbols) both as a function of the parameter eA^. 
Reserving later on the discussion about the time-rescaling with 6, let us observe two 
relevant features. The first one is the presence of a change of curvature of both curves 
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on approaching the critical border eA^ = 6Ii^ ^ 1. The second is the scahng of all points 
in the two curves (one for T1/2S, the other for the remnant magnetization M^.) 

As for the rescaling of the time, let us observe that, due to the particular quasi 
1-D geometry, and to the low concentration S 1, closest dipoles give the major 
contribution to the energy. For instance, the configuration with all spins aligned along 
the 2;-axis will have an energy, 

^'^eA^ (11) 

where the sum is taken over couples {i,j) of neighbor dipoles. In other words 
E' ~ N/{d)^ ~ N6, where {d) is the average distance between two dipoles. 

On the other hand the Landau-Lifshitz-Gilbert equations of motion are invariant 
under a simultaneous scaling of time and energy r' = t/6 and E' = E6 so that we will 
expect r oc 1/6. This simple relation has been verified considering a system with the 
same aspect ratio e, and the same number of particles A^ (so to have the same value of 
eA^) and changing the concentration S over 3 orders of magnitude. Results are presented 
in Fig. [7^, where T1/2 has been shown vs 6. To guide the eye a dashed line indicating 
the inverse proportionality has been superimposed. As one can see, looking at the last 
point to the right side of Fig. [7^, this relation does not hold true for concentrations 
(5 ~ 1, where the nearest neighbor approximation ffTTj) fails. 

The "scale invariance" is even more evident if the average magnetization is 
considered as a function of the rescaled time t6, shown in Fig. [7)d, for the same cases 
belonging to the straight line shown in Fig. [7^). 

At last, we investigate the system behavior on approaching the large A^ limit. First 
of all let us observe that the scaling variable eA^ = RN/ L is well defined in the large A^ 
limit N, L ^ 00, N/L = const and R fixed|§| 

In order to do that we take into account different systems with fixed concentration 
6 and radius R and increasing lenght L and number of particles A^: both in the 
ferromagnetic phase eA^ < 1 ( Fig. [HK), and in the paramagnetic one eA^ > 1 (Fig. [8b) , 
the average magnetization is independent on the number of particles A^. 

5. Conclusions 

In this paper the microcanonical dynamics of a system of random dipoles, interacting 
with a pure dipole-dipole interaction has been considered. We have shown 
that a dimensional "phase" transition, correspondent to a transition from regular 
(ferromagnetic) to stochastic (paramagnetic) regime occurs, in the microcanonical 
ensemble, for low concentration 6. Such transition is characterized from the dynamical 
point of of view by a different behavior of the fluctuations of the average magnetization, 
and from the physical point of view respectively by a zero remnant magnetization, 
Mr = lim^^oo(^2(T)), and finite decay rates, oc I/T1/2S (Paramagnetic Phase) or zero 



§ We thank the Referee for this remark. 
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t5 

Figure 7. a) Dependence of the average demagnetization time T1/2 as a function of 
concentration S, for systems with e = 0.1 and N — 72. The average has been taken over 
an ensemble of 100 different samples. Initially all samples have all spins aligned along 
the z-axis : S'f (0) = 1, i = 1, . . . N . Dashed line represents T1/2 oc 1/5. b) Average 
magnetization {Mz{t)) as a function of the rescaled time t5 for different concentration 
6 and fixed e = 0.1, and = 72 as in a). 

decay rates and finite remnant magnetization (Ferromagnetic Phase). We showed that 
this dimensional transition occurs when the intensive parameter eN = 1, where e is the 
aspect ratio and is the number of dipoles. For instance in an experimental situation 
if we have a non magnetic substrate with R = 1.6 nm and L = 1.6 yum, with lattice 
size ~ 40.4 nm, we expect a dimensional transition for 5 = 0.15%. We also conjectured 
that in correspondence to this transition, the ground state changes from ferromagnetic 
to antiferromagnetic. 

In the future we would like to investigate dilute dipole systems in the canonical 
ensemble, that is letting the system be in contact with a thermal bath. Our analysis 
in the microcanonical ensemble indicated that the behavior of very dilute dipoles in a 
needle geometry is very similar to a 1-D arrays of dipoles. In the 1-D case dipole 
interaction induces a ferromagnetic ground state, and, due to its anisotropy, to a 
breaking of ergodicity [3]. As shown in Ref. [23], the ergodicity breaking threshold 
can induce very large demagnetization times thus producing ferromagnetic behavior in 
finite samples. Thus, even if one would expect that invariant tori will be destroyed 
under a suitable thermal perturbation, the question on the demagnetization times in 
presence of temperature and on the relevance of the ergodicity breaking is still open. 

The ergodicity breaking found in Ref. [3] considers the total magnetization as an 
order parameter. On the other hand, different order parameters can be defined in dipole 



Dynamics of random dipoles : chaos vs ferromagnetism 



14 





1000 



Figure 8. Average magnetization (Mz) vs the dimensionless time r, for different 
sample lengths L, as indicated in the legend, and different number of spins N, at 
fixed density {N/L = 0.36) for the paramagnetic phase Ne — 2.16 (a), and for the 
ferromagnetic one eN ~ 0.1 (b). In a) is i? = 6, 5 = 0.01, while in b) is i? = 4 
and S = 0.0015625. Initially we choose 5f (0) ^ 1, i = 1, . . . , N . An ensemble of 100 
different configurations has been considered. 



systems, depending on the ground state configuration, for instance an anti-ferromagnetic 
order parameter or a spin glass order parameter. Therefore, it would be interesting to 
investigate the existence of an ergodicity breaking energy threshold with respect to 
different order parameters. 

In conclusion dipole-dipole interacting spin systems offer a realistic playground to 
analyze many properties of magnetic systems which challenge our comprehension. 
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